o 

(N 
O 

Q 



C/5 



o 



X 



Jointly Optimal Channel Pairing and Power 
Allocation for Multi-channel Multi-hop 

Relaying 

Mahdi Hajiaghayi, Min Dongf , and Ben Liang 



Abstract 

We study the problem of channel pairing and power allocation in a multi-channel multi-hop relay 
network to enhance the end-to-end data rate. Both amplify-and-forward and decode-and-forward relaying 
strategies are considered. Given fixed power allocation to the channels, we show that channel pairing over 
multiple hops can be decomposed into independent pairing problems at each relay, and a sorted-SNR 



| channel pairing strategy is sum-rate optimal, where each relay pairs its incoming and outgoing channels 

> 



by their SNR order. For the joint optimization of channel pairing and power allocation under both 



total and individual power constraints, we show that the problem can be decomposed into two separate 

ON ■ 

subproblems solved independently. This separation principle is established by observing the equivalence 
between sorting SNRs and sorting channel gains in the jointly optimal solution. It significantly reduces 
the computational complexity in finding the jointly optimal solution. The solution for optimizing power 
allocation is also provided. Numerical results are provided to demonstrate substantial performance gain 
of the jointly optimal solution over some suboptimal alternatives. It is also observed that more gain is 



obtained from optimal channel pairing than optimal power allocation through judiciously exploiting the 
variation among multiple channels. Impact of the variation of channel gain, the number of channels, 
and the number of hops on the performance gain is also studied through numerical examples. 
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I. Introduction 

The emerging next-generation wireless systems adopt a multi-channel relaying architecture 
for broadband access and coverage improvement H), GO. As opposed to a narrow-band single- 
channel relay, a multi-channel relay has access to multiple channels, e.g., different frequency 
channels or subcarriers in an Orthogonal Frequency Division Multiplexing (OFDM) system. It 
may receive a signal from one channel and transmit a processed version of the signal on a 
different channel. This multi-channel relaying capability can be exploited to choose forwarding 
channel adaptively for the incoming signals, taking advantage of the diverse strength of different 
channels. 

In this work, we address the general problem of channel selection and power allocation 
strategies at multi-channel capable relays to forward data in a multi-hop relaying network. This 
problem involves two issues: 1) channel pairing (CP): the pairing of incoming and outgoing 
channels at each relay; 2) power allocation (PA): the determination of power used to transmit 
signals on these channels. In general, for multi-hop relaying, there is strong correlation between 
CP and PA. Intuitively, to maximize the source-destination performance, the choice of CP at 
each relay would affect the choices of CP at other relays, which further depends on the specific 
PA scheme used. The optimal system performance requires joint consideration of CP and PA. 
Our goal is to maximize the end-to-end data rate in a multi-hop relaying network. 

One may view a CP scheme at each relay as a routing scheme embedded in the network router. 
However, despite bearing some resemblance, the CP problem differs from the conventional multi- 
channel routing problem: For channel pairing, the total cost of two paired incoming and outgoing 
links is not additive as it is typically assumed in the routing case. Furthermore, the cost of each 
link cannot be independently defined in CP. The source-destination achievable data rate is dictated 
by the end-to-end signal-to-noise ratio (SNR), which is a nonlinear function of the channel gain 
and power used on each link. 

A. Contributions 

In this paper, we present a comprehensive solution for jointly optimizing CP and PA to 
maximize the source-destination data rate in a multi-channel multi-hop relay network. The main 
results in our work are summarized as follows: 
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• Given fixed power allocation, the sorted-SNR CP scheme is shown to be optimal in multi- 
hop relaying. Specifically, CP can be separated into individual pairing problems at each 
relay, where the relay matches the incoming channels to the outgoing channels in the order 
of SNRs seen over these channels. 

• The problem of joint CP and PA optimization can be decomposed into two separate problems 
which can be solved independently: first CP optimization, and then PA optimization. The 
decoupling of CP and PA optimization significantly reduces the problem search space and 
reduces the complexity of optimal solution. This separation principle holds for both amplify- 
and-forward (AF) and decode-and-forward (DF) relaying strategies, and for either total or 
individual power constraints imposed on the transmitting nodes. 

• In joint CP and PA optimization, the optimal CP is shown to be decoupled into per relay CP 
The channels at two consecutive hops are optimally paired according to their channel gain 
order, without the need for knowledge of power allocation on each channel. This allows 
simple distributed relay implementation for optimal operation, as well as easily adapting to 
the network topology changes. 

• The solution for PA optimization in a multi-hop setting is proposed for both AF and DF 
relaying. For DF relaying, we develop an algorithm through a dual-decomposition approach, 
where we are able to obtain the semi-closed-form PA expression. In addition to depicting 
power distribution across channels at an individual node, the PA expression allows us to 
characterize the interaction among the nodes for power determination on a multi-hop path. 

The separation of joint CP and PA has been established for dual-hop (i.e., single-relay) DF 
relaying in prior work under total power constraints ft3^_ and individual power constraints [f4j . It is 
somewhat surprising that such separation property is preserved in the general multi-hop relaying. 
In fact, the generalization from the dual-hop case to the multi-hop case is non-trivial. For the 
latter, in addition to being a function of power allocation, the pairing at each relay along the 
hops needs to be optimized jointly, adding an additional dimension for the optimization problem. 
Intuitively, to maximize the source-destination rate, the choice of CP at each relay would affect 
the choices of CP at other relays, which also depend on the specific power allocation scheme 
used. Therefore, it is not apparent that the optimal CP can be decomposed into independent 

'A flaw in the proof was later found in the paper. However, with modification, it can be corrected to show the same result. 
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pairing problems at each relay, or that CP and PA can be separately considered. Besides, the 
two different techniques used in j3]|, JH to show the separation result are complicated. They 
cannot be simply extended and applied to the multi-hop case. Instead, we develop a new approach 
to attack the problem that leads to the separation principle of joint CP and PA for both AF and 
DF relaying in a general multi-hop setting, under either total or individual power constraints. 
Our approach provides a rigorous and direct way in proving the separation result. 

We further provide numerical studies on the performance of jointly optimal CP and PA scheme 
and compare it with those of other alternatives for multi-channel multi-hop relaying. We will 
show that, although both CP and PA improves the performance, the optimal CP is more crucial 
than the optimal PA. In other words, a major portion of the gain comes from the optimal CP. In 
addition, we will see that uniform PA with optimal CP achieves near-optimal performance even 
at moderately high SNR, for AF relaying. This significantly simplifies the PA implementation, 
without the need of centralized channel information for either CP or PA. The gain by the optimal 
CP widens with a higher level of channel gain variation across channels, or a larger number of 
channels, indicating that these factors can be judiciously exploited through CP. The optimal PA, 
on the other hand, is insensitive to these changes. Finally, we will also show that the gain of 
jointly optimal CP and PA becomes more pronounced with an increasing number of hops. 

B. Related Work 

For an OFDM system as a typical example of multi-channel systems, the concept of CP was 
first introduced independently in fl5l and iH for a dual-hop AF relaying system where heuristic 
algorithms for pairing based on the order of channel quality were proposed. For relaying without 
the direct source-destination link available, [J5J used integer programming to find the optimal 
pairing that maximizes the sum SNR. From a system-design perspective, the sorted-SNR CP 
scheme was proposed in [6J and was shown optimal for the noise-free relaying case, under the 
assumption of uniform power allocation. 

These works sparked interests for more research in this area. In the absence of the direct source- 
destination link, for the practical case of noisy-relay, by using the property of L-superadditivity 
of the rate function, the authors of proved that the sorted-SNR CP still remains optimal 
for sum-rate maximization in dual-hop AF relaying OFDM system. Subsequently, it was further 
proved in [8], through a different approach, that the sorted-SNR CP scheme is optimal for both 
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AF and DF relaying in the same setup. When the direct source-destination link is available, 
BH presented two suboptimal CP schemes. For the same setup, a low complexity optimal CP 
scheme was later established in IfTOll for dual-hop AF relaying, and the effect of direct path 
on the optimal pairing was characterized. In addition, it was shown in [ITOll that, under certain 
conditions on relay power amplification, among all possible linear processing at the relay, the 
channel pairing is optimal. 

The related problem of optimal PA for a dual-hop OFDM system was studied by many [fTTTl . 
[fT2l . lfT3l for different relay strategies and power constraints. The problem of jointly optimizing 
CP and PA was studied in a dual-hop OFDM system for AF and DF relaying in |[T4|| and 
[fT5l . respectively, where the direct source-destination link was assumed available. The joint 
optimization problems were formulated as mixed integer programs and solved in the Lagrangian 
dual domain. Exact optimality under arbitrary number of channels was not established. Instead, 
by adopting the time-sharing argument [16J in their systems, the proposed solutions were shown 
to be optimal in the limiting case as the number of channels approaches infinity. 

Without the direct source-destination link, jointly optimizing CP and PA for DF relaying in 
a dual-hop OFDM system was investigated in ^ and flU, where 01 assumed a total power 
constraint shared between the source and the relay, and flU considered individual power con- 
straints separately imposed on the source and the relay. In both cases, two-step separate CP and 
PA schemes were proposed and then proved to achieve the jointly optimal solution. For this 
dual-hop setup, it was shown that the optimal CP scheme is the one that maps the channels 
solely based on their channel gains independent of the optimal PA solution. 

Similar studies on the problem of CP and PA in dual-hop AF relaying or multi-hop relaying 
have been scarce. The authors of iTPTll proposed an adaptive PA algorithm to maximize the end- 
to-end rate under the total power constraint in a multi-hop OFDM relaying system. For a similar 
network with DF relaying, |[T8l studied the problem of joint power and time allocation under 
the long-term total power constraint to maximize the end-to-end rate. Furthermore, in [FT], the 
idea of using CP to further enhance the performance was mentioned in addition to PA. However, 
no claim was provided on the optimality of the pairing scheme under the influence of PA. The 
optimal joint CP and PA solution remained unknown. 
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C. Organization 

The rest of this paper is organized as follows. In Section HH we present the system model and 
joint optimization formulation. In Section |inj given a fixed PA solution, we provide the optimal 
CP scheme based on the sorted SNR for both AF and DF strategies. The joint optimization 
problem of CP and PA is considered in Section [TV] where the separation principle between CP 
and PA optimization is established. The optimal PA solution is then discussed in Section |V] for 
multi-hop relaying under both total and individual power constraints. The numerical study are 
provided in Section [VH and finally we conclude the paper in Section IVIII 

II. System Model and Problem Statement 

We consider an M-hop relay network where a source node communicates with a destination 
node via (M — 1) intermediate relay nodes as illustrated in Fig. [Q For broadband communication 
between the nodes, the frequency bandwidth is split into multiple subbands for data transmission. 
A practical system with such an approach is the OFDM system where the bandwidth is divided 
into N equal-bandwidth channels. We denote by h m>n , for m = 1, • • • , M and n — 1, • • • , N, 
the channel response on channel n over hop m. The additive noise at hop m is modeled as 
an i.i.d. zero mean Gaussian random variable with variance a^. We define a mn = 1 ^ as 
the normalized channel gain against the noise power over channel n of hop m. In the rest of 
presentation, we simply refer to it as channel gain without causing confusion. We make the 
common assumption that the full knowledge of global channel gains is available at a central 
controller, which determines the optimal CP and paS. We further assume that the destination is 
out of the transmission zone of the source, and therefore, there is no direct transmission link. 
For M-hop relaying, a transmission from source to destination occupies M equal time slots, 
one for each hop. In the mth slot, m — 1, • • • , M, the mth node (the source node if m — 1, 
otherwise the (m — l)th relay node) transmits a data block to the (m + l)th node (the destination 
node if m = M, otherwise the mth relay node) on each channel. Our study is constrained to 
half-duplex transmissions, where the relay nodes cannot send and receive at the same time on 

2 However, we show later that, for joint CP and PA optimization, the CP solution requires only local channel information at 
each relay, and given the proposed CP solution, a uniform PA scheme without using channel information is near optimal even 
at moderately high SNR for AF relaying. 
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the same frequency. However, the transmission of different data blocks in different hops may 
occur concurrently, depending on the scheduling pattern for spatial reuse of spectrum. 



hop 1 




hop 2 



hop M 




source 



relay 1 



relay (M-l) destination 



Fig. 1. Illustration of multi-channel multi-hop relaying network with channel pairing. 



A. Relaying Strategies 

We consider two types of relaying strategies: AF and DF. In AF relaying, a relay amplifies 
the data received from an incoming channel and directly forwards it to the next node over an 
outgoing channel. In DF relaying, a relay attempts to decode the received data from the previous 
node over each incoming channel and forwards a version of the decoded data on an outgoing 
channel to the next node. We consider the simple repetition-coding based DF relaying |[T9l , 11201 , 
where the relay is required to fully decode the incoming message, re-encodes it with repetition 
coding, and forwards it to the intended user. 

B. Channel Pairing 

The relay conducts CP, matching each incoming channel with an outgoing channel. As different 
channels exhibit various quality, a judicious CP scheme can potentially lead to significant 
improvement in system spectral efficiency. 

We denote path V% = (c(l, i), ■ ■ ■ , c(M, i)), where c(m, i) specifies the index of the channel 
at hop m that belongs to path V%. For example, V{ = (3,4, 2) indicates that path V% consists of 
the third channel at hop 1, the fourth channel at hop 2, and the second channel at hop 3. Once 
channel pairing is determined at all the relays, the total N disjoint paths Vi, • ■ • ,Vn can be 
identified from the source to the destination. 
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C. Power Allocation 

Denote the power allocated to channel n over hop m by P m>n . The SNR obtained on this 
channel is represented by 7 m ri = P m ,n^m,n- For each path V{, let 7 m i = j m ,c( m ,i) represent the 
SNR seen over hop m on this path. 

Let Pj = (Pi )C n,i), • • • , Pm,c(m,i)) be the PA vector for all channels along path V%. The source- 
destination equivalent SNR of path V% is denoted by jsr>(Vi, Pi). For AF relaying, it is given 
by ED, 

^.)=(n(i + ^)-i) 1 . 

and, in Section |Vl we will also use its upper bound ETI . 

whose approximation gap vanishes as the SNR becomes large. For DF relaying, we have 

7sD(^i P i)= min „Tm,i • (3) 
m=l,— ,M 

We consider two types of power constraints: 

a) Total power constraint:: The power assignment P m , n , f° r m = 1 ■ ■ ■ M and n = 1 • • • iV, 
must satisfy the following aggregated power constraint 

M N 
m=l n=l 

b) Individual power constraint:: The power assignment P m ,n, for n = 1, ■ ■ • , N, needs to 
satisfy the power constraint of the individual node m, i.e., 

N 

J2 P m,n = Pmt, m = l,---,M, (5) 

n=l 

where P mt denotes the maximum allowable power at node m. 



D. Objective 

Our goal is to design a jointly optimal CP and PA strategy to maximize the source-destination 
rate under multi-hop relaying. The source-destination rate achieved through path Vi is given by 

RsD(V i ,P l ) = — log 2 (l + 7 SD (7' i ,P < )) J 
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where F s is the spatial reuse factor. In multi-hop relaying that allows concurrent transmissions, 
F s takes value between 2 and M (since F s > 2 under the half-duplex assumption). The sum 
rate of all paths determines the total source-destination rate of the system, denoted as R t , i.e., 



N 



R t = Y, R ™( V ^ p i)- ( 6 ) 



i=l 

It is a function of both {Pi} and {Pi}, which should be jointly optimized: 

max R t (7) 
{ViUPi} 

subject to © or ©, 

PihO, i = l,---,N, (8) 
where >z signifies element-wise inequality. 

III. Optimal Multi-hop channel Pairing Under Fixed Power Allocation 

In this section, we first consider the case when PA is fixed and given. In this case, the 
optimization problem in © can be re-written as 

N 



max 

{■Pi. 



axJ2 R SB(Vi,Pi), (9) 

1=1 

and the optimal CP {V*} is a function of {P^}. To simplify the notation, in this section we 
rewrite Rsui'Pi) an d IsD^i) an d drop their dependency on Pj with the understanding that {P,} 
is fixed. In the following, we solve © to obtain the optimal CP scheme under this fixed PA. 
We emphasize that here the generalization from the dual-hop case to the multi-hop case is non- 
trivial. Intuitively, there is no obvious way to decouple the sequence of pairings at all (M — 1) 
relays. Indeed, the equivalent incoming channel from a source to a relay and the equivalent 
outgoing channel from that relay to the destination depend on how the channels are paired over 
multiple hops. However, we will show that the optimal CP solution over multiple hops can in 
fact be decomposed into (M — 1) independent CP problems, where the mapping of incoming 
and outgoing channels at each relay is only based on the sorted SNR over those channels, and 
therefore can be performed individually per hop. 

In the following, we first establish the optimality of the sorted-SNR CP scheme for the case 
of M = 3 and N = 2, and then we extend the result to arbitrary M and N. 
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source 112 relay 1 122 relay 2 132 destination 



(2) 




Fig. 2. Three-hop relay with two channels. 

A. Optimal Channel Pairing for Three-Hop Relaying 

1) Two-channel case (N = 2): We first consider a three-hop relaying network with two 
channels, as depicted in Fig. [2] Without loss of generality, we assume channel 1 exhibits equal 
or larger SNR than channel 2 over all the three hops, i.e., 



The optimal CP scheme for this case is given in Proposition [IJ 

Proposition 1: For M = 3 and N = 2, the solution to © is the sorted-SNR CP scheme 
performed on each relay, i.e., {V*} = {(1, 1, 1), (2, 2, 2)} under condition (flOl) . 

Proof: The proof essentially exams possible path selections and shows that the sorted-SNR 
per relay provides the highest source-destination sum rate for both AF and DF relaying. See 
Appendix |A] for details. ■ 

2) Multi-channel case (N > 2): Here, we provide an argument to extend the result in 
Proposition Q] to a system with an arbitrary number of channels. 

Proposition 2: For M = 3 and N > 2, the solution to © is the sorted-SNR CP scheme 
performed on each relay. 

Proof: Suppose the optimal pairing does not follow the pairing rule of sorted SNR. There 
is at least one relay (say, Relay 2) that has two pairs of incoming and outgoing channels that are 
mis-matched according to their SNR. That is, there exist two channels i\ and %2 over hop 2, and 
two channels ji and j 2 over hop 3 that are respectively paired with each other while j 2 ,h < l2,i 2 
and 73 > 73 J2 - Note that these two channel pairs belong to two disjoint source-destination 
paths that can be regarded as a 2-channel relay system. From Proposition [fl we know that pairing 
channels i x with j 2 and i 2 with j\ at relay 2 achieves a higher rate than the existing pairing 



7m,l > 7m,2 



for m = 1, 2, 3. 



(10) 
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over these two paths. Hence, by switching to this new pairing while keeping the other paths the 
same, we could increase the total rate. This contradicts our assumption on the optimality of a 
non-sorted SNR CP scheme. Hence, there is no better scheme than sorted-SNR CP to obtain the 
maximum sum rate. ■ 

B. Optimal Channel Pairing for Multi-hop Relaying 

Building on Proposition [2l we next extend the result for 3-hop relaying to a relaying network 
with an arbitrary number of hops (M > 3) in the following proposition. 

Proposition 3: The solution to (J9]) is the sorted-SNR CP scheme individually performed at 
each relay. 

Proof: We prove by induction. It is shown in Proposition |2] that the sorted-SNR CP is 
optimal for M — 3. Suppose the claim holds for M < L. Now consider M = L+ 1 as shown in 
Fig. Oa). Let j eq >n be the received SNR from the source to relay L — 1 over the nth incoming 
channel of that relay. We establish TV equivalent channels between the source and relay L — 1, 
with SNR over the nth channel as 7 e9 , n . Then, the (L+ l)-hop relaying network can be converted 
to a 3-hop network, with an equivalent relay whose incoming channels have SNR {^ eq ,n] and 
outgoing channels remain the same as those of relay L — 1, as shown in Fig. Ob). Hence, from 
Proposition [2l the optimal CP is the one where {^ eq ,n] and {^L,n} are sorted and paired at 
this equivalent relay, and {7l,„} and {7L+i, n } are sorted and paired at relay L. Note that the 
sorted-SNR pairing at relay L is independent of how the channels are paired at the other relays. 

Next, ignore relay L and replace it by equivalent channels from relay L — 1 to the destination. 
We now have a L-hop network. From the induction hypothesis, the sorted-SNR CP is optimal. 
In particular, the incoming and outgoing channels at each of relays 1, 2, . . . , L — 2 are sorted by 
their SNR and paired. Since the SNRs {j e q,n} at the equivalent relay are computed by applying 
©or © over these sorted and paired channels from the source to relay L — l, it is not difficult 
to see that {7i_i, n } and {'jeg^n} are ordered in the same way. Therefore, sorting and pairing 
{leq,n} and {jL,n} at the equivalent relay is the same as sorting and pairing {jL-i,n} and {jL,n} 
at relay L—l. Thus, we conclude that at each of relay 1, • • • , L, the incoming and outgoing 
channels are sorted and paired in order of their SNR. ■ 

The significance of Proposition [3] is that the optimal CP for M-hop relaying is decoupled into 
(M — 1) individual pairing schemes at each relay, each solely based on the SNR of incoming and 
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outgoing channels. This decoupling not only reduces the pairing complexity, but also reveals the 
distributed nature of optimal CP among multiple relays, thus allowing simple implementation 
that can easily adapt to network topology changes. 

Remark: We point out that the existing result of optimal CP strategy for dual-hop relaying is not 
sufficient for the induction to prove Proposition [3l Notice that, in the proof, an M-hop network 
(M > 3) was transformed into an equivalent 3-hop network. Reducing a 3-hop network to a 
dual-hop network would require combining relay nodes with either the source or the destination 
to form an equivalent node and equivalent channel gain. The dual-hop result can only be applied 
to pairing with the equivalent channels, but is not sufficient to show the actual physical channel 
should follow the same pairing strategy. Therefore, Proposition [2] is necessary as the basis to 
prove the general M-hop case. 

In addition, in 0, the L-superadditivity property ||22l is used to show that the sorted-SNR CP 
is optimal in dual-hop AF relaying for sum-rate maximization. That is, if the source-destination 
rate over each path can be shown to be L-superadditive, it follows that sorted-SNR pairing is 
optimal. However, L-superadditivity does not hold for the rate function in general multi-hop 
relaying, where the source-destination rate is a higher dimensional function defined on R M with 
respect to 7 l n , • ■ • ,7Af, n > for a given n. Thus, a similar proof for the optimality of sorted-SNR 
in the dual-hop case is not available to the general multi-hop case. 

IV. Jointly Optimal Channel Pairing and Power Allocation: A Separation 

Principle 

So far, given a fixed PA scheme, we have found that the optimal CP scheme for Q is SNR 
based, which depends on the transmission power allocated to each channel. We next present the 
solution for © by jointly optimizing CP and PA. 

The apparent coupling of CP and PA makes a direct exhaustive search for the jointly optimal 
solution prohibitively complex. Instead, we will show that the joint optimization problem can 
be decoupled into two separate CP and PA subproblems. Specifically, we prove that the jointly 
optimal solution is obtained by pairing channels based on the order of their channel gains 
(normalized against the noise power), followed by optimal PA over the paired channels. This 
separation principle holds for a variety of scenarios, including AF and DF relaying under either 
total or individual power constraints. 
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Fig. 3. Converting an (L+l)-hop relaying to an equivalent 3-hop relyaing. (a) An (L+l)-hop relaying network; (b) An equivalent 
3-hop relaying network. 



Our argument for the separation principle is briefly summarized as follows. We first show that, 
at a global optimum, the channel with a higher channel gain exhibits a larger SNR. This relation 
reveals that the SNR-based ordering of channels is the same as the one based on channel gain. 
Hence, we conclude that the sorted CP scheme based on channel gain is optimal when PA is 
also optimized. 

A. Ordering Equivalence at Optimality 

Let 7j^ n be the received SNR under the optimal PA solution {P*} for hop m and channel 
n. For both total and individual power constraints, the following proposition establishes the 
equivalence between channel-gain ordering and SNR-based ordering at the optimality. 

Proposition 4: In the optimal CP and PA solution for ©, at each hop, the channel with better 
channel gain also provides a higher received SNR, i.e., a mji > a m j implies 7^ > 7^ -, for 
m — 1, ■ ■ ■ , M; i,j G {1, • • • , N}, and i ^ j. 

Proof: We first provide a proof for a multi-hop system consisting of two channels. We then 
explain how it can be extended to a system with an arbitrary number of channels. 

a) N = 2: We prove the proposition by contradiction. Let V\ and V2 represent the two 
disjoint source-destination paths corresponding to the optimal CP scheme. Consider any hop m 
along these paths. Without loss of generality, let channel 1 belong to V2, channel 2 belong to 
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Pi, and a m ,i > a m , 2 . Suppose at optimality 7^ < 7^, i.e., P„ )2 a mi2 > P^ yl a m ^, where P^ 
and 2 are the power allocated to channels 1 and 2, respectively. Let P mt = P^ ll + P^ 2 - 
Consider the following alternate allocation of power between channels 1 and 2 over hop m 

p _ Q m,2 p » p _Sl p * nn 

r m,l — r m,2l r m,2 — r m,l- \ Ll ) 

0>m,\ a m,2 

We further swap the two channels so that channel 1 belongs to path Pi and channel 2 belongs 
to path P 2 . Since P m ,ia mi i = P„. 2 a mi2 and P m|2 a m)2 = Pmi a m,\, the above procedure of power 
re-allocation and channel swapping does not change the end-to-end rate. 

p p _ a m,2 , p p * s Q-m,l p * 

r m,l T f m ,2 — V/mi ■ r m,l/ m,l 

_ Qm,2 p (fl-m,l) 2 — («m,2) 2 p * 

0"m,l Q"m,lU"m,2 

«m,2 p (« m ,i) 2 — (a m ,2) 2 flm,2 p n0 s 

< r mi H j Pnt V iz J 

GWl a mA a m,2 a m,l + ^m,2 



P 



mt 1 



where inequality (fT2l) is obtained from our assumption that P^ x a m ^ < P^ 2 a m2 , which can be 
rewritten as P^ x < P mt , and that a m>i > a TO>2 . This contradicts our initial assumption 

that the original PA is globally optimal. 

b) N > 2: A similar proof by contradiction as it was used in Section ITlI- A2I for Proposition 
[2] can be applied to generalize the above result to N > 2. For an iV-channel relay system with 
N > 2, suppose the optimal CP scheme follows the pairing rule of the sorted CP based only on 
SNR gain and not channel gain. As a result, there is at least one hop over which, between two 
channels, the channel with better channel gain demonstrates a lower SNR. These two channels 
essentially belong to the two source-destination paths that can be considered as a 2-channel relay 
system. From the above, we know that by just swapping these two channels and applying the 
alternate allocation of power in (fTTT) . the sum power is reduced while maintaining the same rate. 
This leads to a contradiction of our early assumption on the optimality of the sorted CP not 
being conducted based on the channel gain. ■ 

B, Separation Principle 

Proposition 5: The joint optimization of CP and PA in © can be separated into the following 
two steps: 
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1) Obtain the optimal CP {V*}. The optimal CP {V*} is independent of {P*} and is 
performed individually at each relay in the order of sorted channel gain. 

2) Obtain the optimal PA {P*} under the optimal CP {V*}: 

N 

{P*} = argmax V R SD (V* , Pi) subject to © or ©. (13) 

{^} 

Proof: From Proposition [3l with optimal PA {P*}, the sorted-SNR CP gives the optimal 
{V*}. From Proposition HI at optimality, the sorted-SNR CP is equivalent to sorting channel 
gains, which does not require the knowledge of {P*}. The optimal {P.*} then can be obtained 
under the optimal CP, and we have the separation principle. ■ 
Decoupling the CP strategy from PA strategy significantly reduces the problem search space. 
In addition, the optimal CP strategy in the presence of multiple hops is further decoupled into 
independent sorting problems at each hop, which only depends on the channel gain on the 
incoming and outgoing channels. The complexity of the optimal CP strategy for each hop is that 
of sorting channel gain, which is 0(N log N). Therefore, the total complexity of the joint CP 
and PA optimization is O (MN log N) in addition to the complexity of PA optimization. 

V. Optimal Power Allocation for Multi-hop Relaying 

So far we have obtained the optimal CP at all relays. We next find the optimal PA solution for 
a given CP scheme as in (fl"3l) . With the channels paired at each relay, the system can be viewed 
as a regular multi-hop system. Without loss of generality, we assume the channel gains at each 
hop are in descending order according to their channel index, i.e., a m> i > a m 2 > • • • > a m ,N, 
for m = 1, • • • , M. From Proposition [5l the channels with the same index are paired, and 
a path with the optimal CP consists of all the same channel index, i.e., V* = («,••• , i). In 
the following, we consider the PA optimization problem for total power and individual power 
constraints separately. 

A. Total Power Constraint 

The optimal PA solution with a total power constraint for a multi-hop relaying OFDM system 
was obtained in ifTTl . The results can be directly applied here. We briefly state the solution for 
completeness. 
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The PA optimization problem in (1131) with a total power constraint has the classical water- 
filling solution 

n 1 1 + 

for % = !,-■■ ,N, (14) 



where [x] + = max(x, 0). The Lagrange multiplier A is chosen such that the power constraint in 
([4]) is met, and a eq ^ is an equivalent channel gain over the path V*. 

For DF relaying, the equivalent channel gain, denoted as a£^, is given by ifTTl 




1 



I > . — I • '' --V. (15) 

In other words, the equivalent channel gain is iV times the harmonic mean of the channel gain 
over each hop. It is obtained following the fact that, to maximize the source-destination rate on 
one path, the total power allocated to the path must be shared among the channels on this path 
such that all channels exhibit the same SNR. The power allocated to each transmitting node on 
path V* is given by 

p* 

a m,i 2-/m'=l ~a~~T- 

For AF relaying, the exact expression for equivalent channel gain on path V* is difficult to 
obtain. However, its upper bound approximation can be expressed as IfTTl 



M - - 2 



a AF « 



' * = V"-,iV. (17) 



\m=l 



In this case, the equivalent channel amplitude (normalized against noise standard deviation) is 
N times the harmonic mean of the channel amplitude over each hop. It is obtained using the 
upper bound approximation of equivalent SNR in © over a path. The power allocated to each 
transmitting node on path V* is given by 

p* 

P m> i = t-^M i (18) 

V a m,i Z^m'=l ran- 

The PA solution in (fT4l) requires global channel gain information and therefore needs to be 
implemented in a centralized fashion. 
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B. Individual Power Constraint 

For DF relaying, the source-destination sum rate in © reduces to 

1 N 

R^ F = ~et min log 2 (l + P m , n a m>n ). (19) 

r s ^— ' m=l,— ,M 
n=l 

Maximizing (fT9l) over {P m>n } under individual power constraints in © can be cast into the 
following optimization problem using a set of auxiliary variables r = [n, • • • ,r7v] T : 

1 * 

max— V] r n (20) 
r,P F s ^ 

n=l 

subject to i) r n < log 2 (l + P m , n a min ), m = 1, • • • , M, n = 1, • ■ ■ , N; 

N 

ii) y^ y P m , n < Pmt, m=l,---,M; 

n=l 

m) P h 0, 

where P = [P m , n ]MxiV- Since the objective function is linear, and all the constraints are convex, 
the optimization problem in (|20l) is convex. For such a problem, Slater's condition holds 11231 , 
and the duality gap is zero. Thus, (1201) can be solved in the Lagrangian dual domain. Since the 
spatial reuse factor F s is a constant, we drop it for simplicity without affecting the optimization 
problem. Consider the Lagrange function for (1201) . 



N N M M / N 

- log 2 (l + P m 

,n(l m.ii)) ^ ^ A m | ^ ^ Pm,n Pmt 
n=l n=l m=l m=l \n=l 

(21) 

where /x = [fJ, m>n ]MxN with |U mjn being the Lagrange multiplier corresponding to constraint (i) 
in (l20l) . and A = [Ai, • • • , Xm] t with A m being the Lagrange multiplier associated with power 
constraint in (in) in (l20l) . The dual function is given by 

g{\ M) = max £(P, r, jx, A) (22) 

r,P 

subject to P y 0. 
Optimizing (|2T|) over r for given P, fi and A yields 

M 

^/i m , n = l, for n=l,---,JV. (23) 



m=l 
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Substituting this into £(P, r, /x, A), we obtain 

N M M 
£(P, r, /x, A) = ^ ^2 (/W lo g2(! + p m ,n^m,n) ^mPm,n) ~\~ ^ ^ ^m-Pmt- (24) 

n=l m=l m=l 

It is clear that the dual function g(n, A) obtained by maximizing (l24l) can be decomposed 
into NM subproblems 

N M M 



n=l m=l m=l 



with 



p* 

m,n 



9mn\f^mm ^m) m&X £ mn (P mn , fJ"m,ni ^m) (25) 

subject to P mj „ > 

where 

for m = 1, • • • , M; n = 1, • • • , N. By applying KKT conditions ll23l to (1251) . the optimal power 
allocation P^ n , as a function of /i m ,n and A m , is derived as 

-r , (26) 

A m Q>m,n _ 

for n = 1, • • • , TV and m = 1, • • • , M, where A m is chosen to meet the power constraint in (|20l) . 
Finally, the optimization problem in (|20l ) is equivalent to the dual problem 

min#(^,A) (27) 

subject to [i >: 0, A y 0; 

^ AW,n = lj for n = l,---,iV. 

m=l 

This dual problem can be efficiently solved by using the projected subgradient method [|24l . 
Analogous to a common subgradient method, a sequence of Lagrange multipliers is generated 
which converges to the optimal A* and /x* minimizing g(/x, A). This convergence is achieved 
provided that a suitable step size is chosen at each iteration [|24l . The difference between projected 
and normal subgradient methods lies in having an extra constraint J2m=i ^m,n = 1- To satisfy 
this constraint, at each iteration, the projected subgradient method projects the columns of /x 
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(obtained by subgradient method) onto a unit space to attain a set of feasible multipliers. At 
each iteration, a subgradient of g(fi, A) at the current values of [i m ,n an d A m is required. Let 
[0», 0\} T denote the subgradient, where M = [9^,- ■ ■ , # MM JV ] T and X = [9 Xl , ■ ■ ■ , 0\ M ] T ■ It 
is obtained from (1241) as 



9nm,n = ^2 (1 + ^m,„am,n) 



(28) 



for m = 1, • • • , M and n = 1, • • • ,N, and 



N 



/ / m,rt' 



71=1 

for m = 1, • • • , M, where P^ n is obtained from (T26l) . 

For completeness, we summarize the projected subgradient algorithm for solving the dual 
problem: 

1) Initialize A (0) and ^°). 

2) Given Am and /iS, n , obtain the optimal values of P^ n in d26l) for all m and n. 



3) Update A (0 through 



A m +1) = Ag-^i/J 



(0 



for m = 1, • • • , M. Similarly, update ^ n followed by unitary space projection 



M r (Z+l) ' 



EM - I, 



where 



rXl+l) = r (i) z/ (/)i + 

P'rrijn L"m-,n u fm,n /i J ' 



(29) 



(30) 



for m = 1, • ■ ■ , M; and n — 1, • • • , iV. i/® and i/jt are the step sizes at the Ith iteration 
for multipliers fi and A, respectively. 
4) Let 1 = 1 + 1; repeat from Step 2 until convergence. 
With the optimal A* and fi*, the optimal power solution P* is determined as in (l26l) . i.e., for 



m = 1, • • ■ ,M and n = 1, 



ClA m Q"m,n 



(3D 



where /i TO „ satisfies the constraint (|23l) , and at the same time, \i* m n and A^ are chosen so that 
the individual power constraints in constraint (ii) of (1201) are met. 
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The expression of P^ n in (1311) provides some insight on the structure of the optimal PA for 
multi-hop DF relaying: For a given fi, the power allocation across channels at each node is 
individually determined following a scaled version of the water-filling approach based on the 
channel gain. The scales are determined jointly among different hops to satisfy the condition 
of /x in (1231) . It essentially requires the received SNR 7 m n at each hop of the same path to be 
equal. 

We now consider the PA problem for AF relaying. Unlike DF, the achievable source-destination 
sum rate for AF is not generally concave in {P m>n }. Therefore, we have a non-convex optimiza- 
tion problem formulated as 

N 

v / 



71=1 

N 



m=l 



m,n lA 'm,n 



(32) 



subject to i) ^2 p ™,n - p ™t, m = 1, • • • , M 



n=l 

a) p y 0. 

To find the PA solution we resort to an upper bound of the sum rate in (|32|) . Based on ©, an 
upper-bound approximation is given by 

i^iflogii+ff;^— ) \ (33) 

n=l \ \m=l ' ' / / 

This upper bound becomes tight as the received SNR P m ,nam,n over each channel increases. 

Lemma 1: R" p in (|33l) is concave with respect to {P m>n }- 

Proof: The proof follows from the concavity of © with respect to {P rn ,n}, which can be 
shown by considering its Hessian matrix. The details are given in Appendix [Bj ■ 

Given Lemma Q3 the optimization of {P m ,n} to maximize R" p is a convex optimization 
problem, and we can solve it in the Lagrangian dual domain using KKT conditions ||23~1 . Although 
a closed-form or semi-closed-form solution for {P rrhn } is difficult to obtain in this case, we can 
solve it numerically via standard convex optimization tools. 

VI. Numerical Results 

In this section, we provide simulation examples to evaluate and compare the performance of 
the optimal joint CP and PA scheme with that of suboptimal CP and PA alternatives. We study 
different factors that affect the performance gap under these schemes. 
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Besides the jointly optimized CP and PA scheme, the following suboptimal schemes are used 
for comparison: 1) Uniform PA with CP: the optimal sorted channel gain based CP is first 
performed. At each transmitting node, the power is uniformly allocated on each subcarrier. 
In addition, for total power constraint, the total power is also uniformly allocated to each 
transmitting node. Therefore, for individual power constraint, P m ,n = an d for total power 
constraint, P m>n = -p-^; 2) Opt. PA without CP: only power allocation is optimized but no 
pairing, i.e., the same incoming and outgoing channels are assumed; 3) Uniform PA without CP: 
the same incoming and outgoing channels are assumed, then uniform PA as in the case 1 is 
used. 

We use an OFDM system as an example of a multi-channel system, and refer each subcarrier 
as a channel in this case. For the multi-hop setup, equal distance is assumed from hop to hop, 
and is denoted by d r . No direct link between source and destination is available. The spatial 
reuse factor is set to F s = 3 (i.e., interference is assumed negligible three hops away). We 
assume M = 4, unless it is otherwise specified. An L-tap frequency-selective fading channel 
is assumed for each hop. We define the average SNR as the average received SNR over each 
subcarrier at each receiving node under uniform power allocation. Specifically, it is defined for 
different power constraint as follows: under the total power constraint, SNR avg = £^2, where a 



denotes the pathloss exponent and a 2 the noise variance; under the individual power constraint, 

ml u>r 

Na 2 ■ 



SNR aw - 



A. Impact of the average SNR 

We compare the the performance of various CP and PA schemes at different average SNR 
levels. Fig. H] shows the normalized source-to-destination per- subcarrier rate vs. the average SNR, 
for DF relaying under the total power constraint. The number of channels is set to iV = 64. 
We observe that joint optimization of CP and PA provides significant performance improvement 
over the other schemes. In particular, compared with uniform PA without CP, the optimal CP 
alone provides 4dB gain, and subsequently optimally allocating power provides an additional 
1.5-2dB gain. Interestingly, it is evident that channel pairing alone provides more performance 
gain than power allocation alone does. 

Fig. \5\ plots the normalized source-to-destination per-subcarrier rate vs. the average SNR for 
AF relaying under the total power constraint. Again, iV = 64 is used. For schemes with PA 
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optimization, the upper-bound i?" p in (1331) is used to obtain the PA solution. The actual rate R t 
obtained (as in the objective function in (1321) ) with such PA solution provides a lower bound on 
the rate under the optimal PA. In Fig. [5j for the jointly optimal CP and PA scheme and optimal 
PA without CP scheme, we plot both upper bound and lower bound of the rate for the optimal 
PA solution. We see that these two bounds become tighter as the average SNR increases, due 
to the improving accuracy of approximation R" p . The PA solution derived using i?" p becomes 
near optimal. 

Comparing the performance of different CP and PA schemes shown in Fig. [5j it is seen that 
similarly as in the DF case, joint optimization of CP and PA provides noticeable improvement 
over the other schemes. The gain mainly comes from choosing CP optimally, which provides 
around 2dB gain over no CP schemes. We further observe that, with optimal CP, the gap between 
optimal and uniform PA vanishes at higher SNR, indicating that uniform PA achieves the optimal 
performance at a moderately high SNR range (around 15dB). Interestingly, this is not the case 
for the schemes without CP. The intuition behinds this is the following: At relatively high SNR, 
it is known that the water-filling PA solution in (fl4)) approaches a uniform allocation. Thus, 
the total power is approximately equally distributed to different paths. The power on each path 
P* is then further assigned optimally to each channel on the path according to (fl"8l) . which is 
typically not uniform. The exception is when each hop exhibits a similar channel gain. This 
is more likely to occur as a result of channel pairing, where channels with the same rank, 
more likely with similar strength, are paired with each other. Therefore, with CP, the optimal 
PA approaches to a uniform allocation at a faster rate with increasing SNR This interesting 
observation suggests that, because of CP, at moderately high SNR, we are able to reduce the 
centralized PA solution to a simple uniform PA which requires no global channel information 
without losing much optimality. Note that the same argument is applicable to DF relaying, but 
the optimal PA approaches to a uniform allocation at a much slower rate than that for AF, which 
can be shown by comparing (fl6l) and (fl"8l) . The range of SNR values under consideration is too 
small to see the same effect in Fig. |4l 

Under individual power constraints, the performance comparison of CP and PA schemes are 

3 Note that, for water-filling PA, as SNR — > oo, it approaches to a uniform allocation in all schemes with or without CP. The 
difference is the rate at which PA approaches to a uniform allocation. 
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given in Figs. [6] and [7] for DF and AF relaying, respectively. We assume N = 16. These figures 
further demonstrate the significant improvement by jointly optimizing CP and PA, where most 
of the gain comes from optimal CP. In addition, under AF relaying, we again observe a near- 
optimal performance by uniform PA with CP at moderately high SNR. This suggests that, under 
individual power constraints, the optimal PA is close to a uniform allocation at high SNR as 
well when CP is adopted. This potentially simplifies greatly the PA implementation to achieve 
the optimal performance. 

B. Impact of the Variation of Channel Gain 

In this experiment, we show how the level of channel gain variation across N channels 
affects the performance of various CP and PA schemes. Towards this goal, we increase the 
number of taps of the time-domain frequency-selective channel {i.e., the maximum delay of the 
frequency-selective channel). This increases the level of variation of the corresponding frequency 
response. Figs. [8] and [9] plot the normalized per-subcarrier rate vs. the number of taps of the 
frequency-selective channel for DF and AF relaying, respectively. The number of subcarriers 
is set to iV = 64 and SNR avg = 12dB. As we see, the performance gap between the schemes 
with optimal CP and without CP increases as the level of channel gain variation increases. This 
demonstrates that the optimal CP schemes benefit from an increased level of channel diversity, 
which is utilized effectively through the channel pairing. On the other hand, the relative gain of 
optimal PA to uniform PA is insensitive to such change and remains constant. 

C. Impact of the Number of Channels 

In this experiment, we examine the effect of the number of channels, under the same level 
of channel gain variation across channels, on the performance of various CP and PA schemes. 
For different N, the subcarrier spacing (i.e., bandwidth of each channel) is fixed. In order to 
set the same level of channel gain variation in frequency, we keep the maximum delay of the 
time-domain frequency- selective channel unchanged. Figs . [101 and ITTI demonstrate the normalized 
per-subcarrier rate with respect to N for DF and AF relaying, respectively. The average SNR 
is set to SNR avg = 12dB. We observes that the gap between the two sets of schemes, with and 
without CP, widens as the number of channels increases. The reason behind this observation is 
that, as more channels becomes available, they can be exploited more judiciously for pairing, 
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and therefore, more gain is achieved by CP. The different PA schemes are not sensitive to the 
change of N. 

D. Impact of the Number of Hops 

In this experiment, we study how the number of hops affects the performance of various CP 
and PA schemes. For this purpose, we increase the number of hops while keeping the distance 
between each hop unchanged. Figs. [[2] and [13] illustrate the normalized per-subcarrier rate vs. 
the number of hops with total power constraint for DF and AF relaying, respectively. We set 
N = 64 and SNR avg = 12dB. As expected, for all schemes, the normalized per-subcarrier rate 
decreases as the number of hops increases. For DF, this is because on average the minimum rate 
among all hops decreases as the number of hops increases; for AF, the rate decreases due to 
noise amplification over hops. Comparing different schemes, we observe that the performance of 
the jointly optimized CP and PA scheme has the slowest decay rate, and the performance of the 
schemes with CP decay is slower than those without CP. In other words, the gain of optimal CP 
and PA is more pronounced as the number of hops increases. A multiple-fold gain is observed 
at a higher number of hops. 

VII. Conclusion 

In this paper, we have studied the problem of jointly optimizing spectrum and power allocation 
to maximize the source-to-destination sum rate for a multi-channel M-hop relaying network. For 
fixed power allocation, we have shown that the general CP problem over multiple hops can be 
decomposed into (M — 1) independent CP problems at each relay, where the sorted-SNR CP 
scheme is optimal. We then proved that a jointly optimal solution for the CP and PA problems 
can be achieved by decomposing the original problem into two separate CP and PA problems 
solved independently. The solution obtained through the separate optimization bears considerably 
lower computational complexity compared with exhaustive-method alternatives. The separation 
principle was shown to hold for a variety of scenarios including AF and DF relaying strategies 
under either total or individual power constraints. For all these scenarios, the optimal CP scheme 
maps the channels according to their channel gain order, independent of the optimal PA solution. 
Finally, the solution for PA optimization under the individual power constraints is derived for 
both AF and DF relaying. Significant gains in data rate were demonstrated by employing jointly 
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optimal CP and PA in multi-channel multi-hop relaying. It was also observed that more gain 
is obtained from optimal CP than optimal PA through judiciously exploiting variation among 
multiple channels. 

Appendix A 
Proof of Proposition [j] 

At relay 1, there are two ways to pair the channels: (1) channels 1 and 2 over hop 1 are matched 
with channels 1 and 2 over hop 2, respectively; (2) channels 1 and 2 over hop 1 are matched with 
channels 2 and 1 over hop 2, respectively. These two ways of pairing lead to the following two sets 
of disjoint paths from the source to the destination: {V- 1 ^} = {(1, 1, c(3, 1)), (2, 2, c(3, 2))} and 
{Pf ] } = {(l,2,c(3, 1)), (2, l,c(3,2))}, where the superscript j in {p\ j) } indicates a different 
set of path selection. 

By considering the equivalent channels from the source to the second relay, using the existing 
optimality result for dual-hop relaying [8), it is easy to see that c(3, 1) = 1 and c(3, 2) = 2 are 
optimal for {V^ }. Furthermore, we only need to show 

log 2 (l + TsdC^)) + log 2 (l + 7sd (VP)) > 

log 2 (l + 7sd (V{ 2) )) + log 2 (l + 7sd (Vf )) , (34) 

for the case of c(3, 1) = 1 and c(3,2) = 2 for both {P[ • } and {V^}, since the case of 
c(3, 1) = 2 and c(3,2) = 1 for {pi 2) } can be similarly proven. Inequality (|34|) for the AF and 
DF relaying cases are separately proven as follows: 

a) AF Relaying: By inserting (OQ) into inequality (1341) we need to show 

(i + (Q? - i)- 1 ) (i + (Q? - 1)- 1 ) > 

(l + (QS 2) - l)- 1 ) (l + {Qf - l)- 1 ) , (35) 
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where 



(i) 



(2) 



Ql 



l 
l 

i + 
i + 



l 

7i,i 
1 

71,2 

1 

7i,i 
1 

71,2 



1 + 



1 

72,1 

1 

72,2 
1 

72,2 
1 

72,1 



1 + 
1 + 
1 + 
1 + 



1 

73,1 

1 

73,2 

1 
73,1 

1 

73,2 



The following lemma is used to prove (|35l) 
Lemma 2: With condition (flOl ), we have 



(gS 1) -i)(g^-i)<(Qi ; 



(i) 



)(2) 



Proof: By substituting ((361) in the following term and expanding it, we have 



(Q 



Q? + Q? 
l 

i + — 

7i,i 
1 

1 + 

7i,i 

1 1 

72,1 72,2 



l)-(Qi 2) -l)(Qr J -l) 



(2) 



Qi ^ Q2 



1 

72,2 
1 

72,1 

1 + 



1 

1 + — 

73,1 

1 

1 + — 

73,1 



73,2 



1 



1 

7l,2 



1 

7l,2 

1 

7l,2 

1 + 



1 + 
1 + 



1 

73,1 



1 

72,1 

1 

72,2 

1 + 



1 

73,2 
1 

73,2 



1 

7i,i 



< 



(36) 



(37) 



(38) 



(39) 



where we have used the fact that Q± Q% = Q^Q^ t0 arrive at (|38l) . From condition (fTOb , the 
first product term in (1391) is negative and the second product term is positive, and therefore we 
obtain the last inequality. ■ 
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Consider the subtraction of the RHS from the LHS of (1351) . 
LHS of ([35]) - RHS of d35]) 

= ((q? - 1)- 1 + (q? } - 1)- 1 + (q? - - 1; 



-i 



,4 

-1 



((Ql 2) - ir 1 + (Q? } - ir 1 + - ir\Q? - i 



> A{Q { ? - 1){Q { ? - I) -BiQ?- 1)^-1) (40) 
>to , n( 2 ) n( 2 ) 



1 \ / 1 \ / 1 \ / 1 
1 + — 



— Q2 + Qi Qi Q2 

\72,2 72,1/ LV 71,2/ V 73,2/ V 7l,l/ V 73,1, 

> 0, (41) 

where the inequality (|40l ) holds because of Lemma [2} and the fact that Q^' — 1 > 0, for i — 1, 2 
and j = 1, 2, 3; and the inequality (|4TT) holds because of condition (TTOT ). 
Z?) Di 7 Relaying: Inserting ([3]) into inequality (|34|) . we need to show 

(1 + min(7i 5 i, 72,1, 73,1)) (1 + min(7i )2 , 72,2, 73,2)) > 

(1 + min(7i ) i, 72,2, 73,1)) (1 + nxin(7i )2 , 72,1, 73,2)) • (42) 

We can verify (|42l) by enumerating all possible relations among 7 min , for all m = 1,2,3 and 
n = 1,2, subject to condition (flOl) . For example, when 7^ < 7 2) i < 7 3) i, 7i j2 < 72,2 < 73,2, 
72,2 < 71,1 < 73,2, and 73,2 < 72,1, (52]) reduces to 

(1 + 7i,0(l + 71,2) > (1 + 72,2) (1 + 71,2)- 



The above inequality clearly holds based on the assumption of {7ij} relations. Inequality ((42]) 
can be similarly verified for all other {7^} relations. The details are omitted for brevity. 

Appendix B 
Proof of Lemma Q] 

Let R n denote the end-to-end data rate on path n, we have R n = y- log 2 (1 + j n ), n = 
1, • • • ,N, where ^ 

i^lEir 2 -) • («) 



\va=l 



m,nu"m,n 
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Then R" p = J2 n =i ^° snow -R" P is concave in {P m , n }, it suffices to show that each R n is 
concave in {P mn }. The concavity proof of R„ follows the concavity of 7 n due to the composition 
rules which preserve concavity ||23l . For simplicity, we drop the subscript n from notations in 
(|43T) . In the following we prove that j(P) is concave in P, where P = [Pi, ■ ■ ■ , Pm] t - The 
second-order partial derivatives of j(P) are given by 

i)--{I>) ' / W 



and 




dP 3 dP k a 3 Pf a k P t 



(44) 



for k ^ j. (45) 



»=i 



Hence, the Hessian matrix A 2r y(P) can be expressed as 

A2 ^> = it At) ' (-* (t ±) *ag (^■■,-± sr )+**f), < 46 > 



J=l 



a iPi J K^lPf ' a MPli . 

where q = \q 1 , ■ ■ ■ , q M ] T with q m = — W, m = 1, • • • ,M, and diag(x) denotes the diagonal 

0"m,r m v ^ 

matrix with diagonal elements being the elements in vector x. To prove concavity, we need to 
show A 2 7(F) -< 0. For any vector v = [vi, ■ ■ ■ , v M ] T , we have 



M „ \ 3 / / M , \ M / M 




< 



where the inequality is obtained by using the Cauchy-Schwartz inequality (e T e)(c T c) > (e T c) 2 
for two vectors e = [ei, • ■ ■ , euY and c = [ci, ■ ■ ■ , cm] t , with = -, = and q = /' . 

\ " ' y II, I; 

Therefore, A 2 7(F) ^ 0. 
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2.2 




Avg SNR(dB) on each Channel 



Fig. 4. Normalized rate vs. the average SNR for DF OFDM relaying with M — 4 and N — 64 under total power constraint. 




Fig. 5. Normalized rate vs. the average SNR for AF OFDM relaying with M — 4 and N — 64 under total power constraint. 
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Avg SNR(dB) on each Channel 



Fig. 6. Normalized rate vs. the average SNR for DF OFDM relaying with M — 4 and N = 16 under individual power 
constraint. 



-D- Opt joint PA and CP, UB 

- □ - Opt joint PA and CP, LB 

- -"-Uniform PA w/ CP 
-e- Opt PA w/oCP, UB 

- 8 - Opt PA w/o CP, LB 

- — »- Uniform PA w/o CP 


j&S A _ 

/* / SB 
S*Y S> 

s*Y Y' »l 
A Y /* 7* 

Y* Y Y' Y 
Y, Y Y' Y 
Y*Y Y* Y. 
/ Y> Y 




Y' Y Y' Y 

Y* Y' y 
n> /> Y 

■ /** ■ ■ y~ : ~ 


rf/ 


YV s&' Y 

Y^ * sY 



Y 1 1 1 1 1 1 ' 1 1 

4 6 8 10 12 14 16 18 20 22 

Avg SNR(dB) on each Channel 



Fig. 7. Normalized rate vs. the average SNR for AF OFDM relaying with M — 4 and N — 16 under individual power 
constraint. 
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Fig. 8. Normalized rate vs. number of taps for DF OFDM relaying with M = 4, N = 64, and SNR avg = 12dB under total 
power constraint. 
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Fig. 9. Normalized rate vs. number of taps for AF OFDM relaying with M = 4, N = 64, and SNR avg = 12dB under total 
power constraint. 
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Fig. 10. Normalized rate vs. number of channels for DF OFDM relaying with M = 4 and SNR avg = 12dB under total power 
constraint. 
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Fig. 11. Normalized rate vs. number of channels for AF OFDM relaying with M = 4 and SNR avg = 12dB under total power 
constraint. 
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Fig. 12. Normalized rate vs. number of hops for DF OFDM relaying with N = 64 and SNR avg = 12dB under total power 
constraint. 
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Fig. 13. Normalized rate vs. number of hops for AF OFDM relaying with N = 64 and SNR avg = 12dB under total power 
constraint. 
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